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Hadron spectroscopy on dynamical configurations are faced with the difficulties of dealing with 
the mixing of single particle states and multi-hadron states (for large spatial volumes and light 
dynamical quarks masses). It is conceivable that explicit multi-hadron interpolating operators will 
be necessary for obtaining sufficiently good overlap on to multi-hadron states in order to extract 
the low-lying excitation spectrum. We explore here the feasibility of using four noise diluted 
all-to-all quark propagators in the construction of explicit two-hadron operators on quenched, 
anisotropic lattices. Our longer term goal is to use these operators on large anisotropic, dynamical 
configurations for hadron spectroscopy. 
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1. Introduction 

The goal of the Hadron Spectrum Collaboration is to determine the low-lying hadron spectrum 
from first principles lattice QCD simulations. The low-lying spectra of baryons in the quenched 
approximation has been reported in Refs. [|J]. There has been progress made in dynamical simula- 
tions in the last few years by many collaborations as more and more light, dynamical confurations 
had become available (See [Q] for recent reviews). 

One of the fundamental challenges in spectroscopy of excited states is the issue of identifying 
multi-particle states in the spectrum. It is our expectation that explicit multi-particle operators 
will become necessary to determine the spectrum above multi-particle thresholds on (large) light, 
dynamical configurations. The shift in energy in a finite volume of two-particle states ^ may 
suggest that finite momentum operators may have better overlap when the interaction between the 
two particles is not negligible. Simulation of correlation matrices with finite momenta operators 
using point-to-all quark propagators do not appear to be more cost effective than using stochastic 
estimates for all-to-all propagators with some improvement scheme. Moreover, the number of 
contractions increase dramatically for two-particle correlation functions constructed from point- 
to-all propagators which may result in a significant fraction of the simulation time being used for 
analysis. For these reasons, it is favourable to use some stochastic estimate of all-to-all quark 



propagators (Q9Q-Q35I]) in the simulation of low-lying multi -particle spectra. 



In this study, we continue our efforts on examing the efficacy of the "dilution method" ([ 1C , 22, 



30]) for estimating all-to-all quark propagators for the purpose of simulating multi-hadron states. 
Earlier works have demonstrated the efficacy of this method for zero momentum meson operators 
and group-theoretically motivated baryon operators [f39|]. A detailed study of the dilution scheme 



dependence of three-quark correlators has been reported in these proceedings [|40j]. We shall focus 
here on finite momentum pion (two-quark) operators and k + k + (four-quark) operators. 



2. Construction of Operators/Correlators 
2.1 Parameters 

We have used 100 quenched, anisotropic {a s /a t = 3) Wilson gauge lattices (12 3 x 48) at j8 = 
6.1 and quark masses which correspond to pion masses of approximately 700 MeV. The gauge 



fields were stout-smeared [38] with (n,e) = (16,0.15625) and the quark fields were further smeared 



using gauge invariant Gaussian smearing with (n,a) = (32,3.0). 

2.2 Two-Quark Operators 

Two independent Z4 noise sources 77^1 (jc,f) and their solutions 0m (-?>*) (with A = 0, 1) were 
generated in order to construct the meson correlation functions. A meson correlation function is 
given by, 

C(p,t)= £ [e- iH "-" o) Q{x,t + t ;x ,t pQ{x ^x,t + t )r 
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Figure 1: The pion effective mass using time-spin- 
colour diluted quark propagators. A single time- 
slice was used for the source operator and a peri- 
odic definition of the effective mass was used with 
Af = 3. 
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Figure 2: The pion effective masses for the first 6 
lowest momenta with time-spin-space-diluted quark 
propagators. A single time-slice was used for the 
source operator and a periodic definition of the ef- 
fective mass was used with At = 3. 



where Q(x,t + ?o;*o,^o) is the quark propagator from (xo,to) to (x,t + to). In terms of the noise 
sources and solutions, this correlation function becomes 



C(p,t) 



I 

to,xo,x 



-ip-(x-xo) 



(y o] r</>[i]) (2,0 ( T 7[ t i] rt 0[o]) (*o,to) 



where the meson creation and annihilation operators (in parentheses) are contracted over the colour 
indices. This construction requires Dirac matrix inversions for each noise index A since Tj[ A ] (x,t) 
must be given on all time-slices (and their corresponding solutions). This would have had the ben- 
efit of increased statistics as one has N t measurements on each lattice, but can be computationally 
demanding because the number of inversions increases by a factor of 48 (for this lattice). How- 
ever, one can settle with the source on a single time-slice only, to, if we use 75 Hermiticity for the 
anti-quark propagator. In other words, we can write 



C(p,t) 



to,*o,x 



0[o]^ r 0[i] ) (2,0 (n[i](fciy*?[o] ) (2o,?o) 



for the meson correlation function, which only contains the random noise source at to. The pion 
effective mass using time-spin-colour dilution with this method is shown in Fig. [l[ The cost of 
producing such an effective mass is twice as much as that for traditional point-to-all propagators 
since two independent noise sources must be used in the stochastic case. 

It is just as simple to form correlation functions with finite momenta as with zero momentum 
since the correlation function is a product of meson creation and annihilation operators which can 
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have any of the allowed momenta values. We construct the finite momenta operators given by, 

X 
X 

from which we make the finite p correlation functions (post-hoc). The effective pion mass for 
the six lowest momenta are shown in Fig. 22 where spin+space-cubic dilution was used for the 
stochastic estimates of the quark propagators. It is clear that one can measure up to (at least) 
five distinct momenta with 100 configurations. Note that the cost comparison with the point-to-all 
method is difficult here as this depends on how many source vectors were chosen for the creation 
operator in the point-to-all simulation. 

2.3 Three-Quark Operators 

The dilution test for three-quark operators have been reported in [[39]] and p0|]. Diagonalization 



with several different nucleon-type operators was found to be stable at relatively low dilution levels. 

2.4 Four-Quark Operators 

We have computed the I = 2,%% correlation function as an example of an operator involving 
four-quarks. The correlation function for this state has the form, 

C;nr(f,fo) = ^[0,1] ^[2,3] ~~ ^[0,1,2,3] 

where Ct^L (Mo) is the pion correlation function constructed from noise sources A and B and the 
"crossed" diagram is given by 

C [0 ,i,2,3] {to + t,t Q )= M^ >1} (t )Mf; k 2] (t + t)M^ 3] (t )M l £ 0] (t + 1) 
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We have used 'zero-momentum' pion operators in this study, but any of the other finite momenta 
operators could have been used as well. (Note that the phase shift could also had been determined 
by diagonalizing the matrix of correlation functions using these 'finite momenta' operators.) An 
example of the four-quark correlation function with spin and space-cubic dilution is shown in Fig. ||. 

3. Simulation/Results 

3.1 Parameters 

For this exploratory study, we have used the same quenched (100) configurations that were 
used in the excited baryon spectrum study [fin. These are the quenched, anisotropic (a s /a t = 3) 
Wilson action lattices at /3 = 6.1, 12 3 x 48. The spatial lattice spacing is 0.1 fm. Anisotropic 
Wilson fermions with pions masses of roughly 700 MeV were used for the valence quarks. 

3.2 Two-Quark Operators 

The dilution scheme dependence of the error of the pion mass is shown in Fig. ||. The error 
is determined from a fit to the pion correlation function starting at t m i„ = 10. This value of t m i„ 
was not necessarily the optimal choice for all of the dilution schemes (in terms of % 2 etc) but a 
common choice is shown here for a lucid comparison. The dashed line is the line 1 / y/N, which 
is the expected behaviour for N repeated measurements with N different noise pairs. The fig- 
ure indicates that there may be better dilution choices than the ones investigated here (time+spin, 
time+spin+colour, time+spin+space-cubic). 

A similar plot for pions with finite momentum (p = (1,1,0)) is shown in Fig. |6| The effect of 
colour dilution and in particular, the spatial-cubic dilution is self-evident. This result may have been 
anticipated, since the projection onto this particular momenta is not contaminated by the random 
noise when spatial (cubic) dilution is used. 

3.3 Four-Quark Operators 

The dilution scheme dependence of the error of the n + 7i + correlation function is shown in 
Fig. Q. The error in this case was the average of the errors in the region of t where plateau for 
various dilution schemes were observed/expected. (We have used 5 < t < 19 in this example.) One 
can see from Fig. ^| that space-dilution appears to be an effective dilution for scattering states, just 
as it was for pions with finite momenta. The cost comparison to the point-to-all case here is also 
difficult as the four-quark contraction is nontrivial there, but it is a factor of 4 more as far as the 
inversions are concerned for the stochastic case since there are four independent noise sources. 

4. Summary 

The dilution method for approximating all-to-all quark propagators has been explored for 
two-quark operators with finite momentum and the simplest four-quark state (7i + n + ). We were 
able to obtain a signal for the six lowest momenta states for the pion with 100 configurations and 
spin+space-cubic dilution. Using the same operators that were used to measure the pion correla- 
tion function (i.e. without further simulations), we have demonstrated that a clean signal for the 
simplest four-quark state can be obtained with spin+space-cubic dilution. 
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Figure 5: The error of the pion mass as a function 
of 1 /\/N where N is the number of quark inversions 
needed in each dilution scheme. 
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Figure 6: The error of the pion energy with mo- 
menta (1,1,0) (in lattice units) as a function of 
1 / V^V where N is the number of quark inversions 
needed in each dilution scheme. 



Simulations for five-quark operators are underway as well as the dynamical simulations of 
hadronic operators involving up to five quarks. A new all-to-all method for quark propagators 
which do not rely on stochastic estimates is also being explored. 
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